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Abstract ■ The APT system, its stage of develop- 
ment, hardware, and operations are described. The 
algorithms required to perform the real-time func- 
tions of navigation and profiling are presented. 

The results of computer simulations demonstrate the 
feasibility of APT for its primary mission: topo- 

graphic mapping with an accuracy of 15 cm in the 
vertical. Also discussed is the suitability of 
modifying APT for the purpose of making vertical 
crustal movement measurements accurate to 2 cm in 
the vertical, and at least marginal feasibility is 
indicated. 

Introduction 

The Charles Stark Draper Laboratory, Inc., under 
sponsorship of the U.S. Department of the Interior, 
Geological Survey*, is developing an airborne- 
instrument system capable of providing precisely 
located, geophysical data. This equipment, called 
APTS (Aerial Profiling of Terrain System), is de- 
signed to be carried by a relatively small maneu- 
verable aircraft. It incorporates a surveying 
instrument package, unaffected by aircraft motion, 
capable of continuously providing a high-accuracy, 
three-coordinate position datum. A self-contained 
terrain-measuring device for recording the vertical 
position of the terrain below the aircraft is also 
provided for applications such as stream-valley 
profiling. A vidicon camera is to be boresighted 
with the terrain-measuring device to aid in the 
data processing. Performance goals call for locat- 
ing an unsurveyed ground control point with respect 
to three or more established ground control points 
to an accuracy of ±15 cm in the vertical coordinate 
and ±60 cm in the horizontal coordinate. 

The work completed to date consists of engineer- 
ing analysis, system configuration design, mathe- 
matical analysis, and computer simulation, selec- 
tion and specification of all major hardware com- 
ponents, mechanical design, thermal design and test, 
electronic design, subsystem interface design and 
specification, selection and recommendation of 
flight computer, analysis of flight strategy, spec- 
ification and design of key flight algorithms, and, 
in addition, a flight test of a breadboard laser 
profiler. Within the time and funding resources 
available, the program has been organized so that a 
balanced effort has been maintained in the various 
aspects of system and subsystem design to facilitate 
an orderly transition into the fabrication, labora- 
tory test, and flight trial of a test bed system. 

The APT system is designed to function in a small 
aircraft at altitudes of 1,000 meters or less. The 
system is to be used for applications such as: 

(1) Producing topographic maps. 

(2) Testing the reliability of older published 
maps. 

(3) Fitting specified flood magnitudes into 
local stream-valley geometry i 

Proc. of the 9th GEOP Conference, An International Symposium on the Applications of 
Geodesy to Geodynamics. October 2-5. 1978. Dept, of Geodetic Science Rept. No. 280, The 
Ohio State Univ.. Columbus, Ohio 43210. 


(4) Classifying public lands for waterpower 
potential. 

(5) Mapping the earth's gravity field. 

(6) Definition of temporal change in groundpoint 
elevation or position: 

(a) Subsidence in and around producing oil 
fields, geothermal reservoirs, and heav- 
ily pumped areas of continuing water 
withdrawal. 

(b) Open-pit and strip-mine mineral extraction 
and land reclamation. 

(c) Volcano inflation, faulting, landslides, 
and beach and slope erosion. 

APT System Configuration 

The airborne- instrument package (Figure 1) con- 
sists of a three-gimbaled inertial platform to de- 
fine the position of the aircraft in three coordi- 
nates. A two-axis laser tracker is mounted on the 
same base as the inertial platform in order to up- 
date the long-term drift of the inertial platform. 
Three or more surveyed retroref lectors interspersed 
with several unlocated retroref lectors provide 
ground truth. The inertial platform and laser 
tracker provide the high-accuracy position datum. 

A laser profiler, to provide line scans of the 
terrain, is provided for those applications re- 
quiring terrain profiling information. However, 
the three-coordinate reference system may be used 
with other systems such as scanning lasers or 
radars, side-looking radar, aerial cameras, and 
radiometers. The laser altimeter or profiler is 
just one of a number of equipments that may be used 
in conjunction with the position and attitude meas- 
urement system. An airborne digital computer ac- 
cepts data from the three sensors and performs the 
necessary computations and data processing for 
alignment of the IMU, and position or velocity up- 
dating calculations from range and angles. In 
addition, the computer outputs data to the 
magnetic-tape recorder and display unit. 

The initial flight test configuration will dif- 
fer somewhat from Figure 1. The purpose of the 
initial flights is to basically prove the concept, 
i.e., locate the position of unsurveyed retro- 
reflectors with respect to three surveyed retro- 
reflectors. The basic differences in this initial 
flight configuration compared to the configuration 
depicted in Figure 1 include the fact that the pro- 
filer will be hard-mounted (the angular errors 
from the profiler laser beam will be computationally 
corrected) . Interfaces for leveling the profiler 
platform and reading its gimbal angles are provided, 
however, so that a stabilized profiler platform may 
be added in the future. Another difference in the 
initial flight configuration is that the autopilot/ 
guidance function indicated in Figure 1 will not be 
implemented. Finally, the indicated Kalman filter 
update, combining tracker and IMU data, will not be 
performed in real time. Tracker and IMU data will 


Under contract 14-08-0001-14578 


171 




be recorded in flight, and filtering of the data 
will be performed post-flight. A simple reset of 
position and velocity will be performed in real 
time, however, so that the acquisition time of a 
new retroref lector is minimized. 

A mechanical schematic (Figure 2) of the system 
shows the geometric relationship and functions of 
the gimbaled subassemblies. The IMU and tracker 
are mounted in a common housing. The profiler- 
vertical camera assembly is hard mounted to the 
aircraft, but referenced to the IMU pitch and roll 
isolation axes to account for motion of the air- 
craft. 

The IMU consists of a stable member which houses 
three high-performance gyros and three acceler- 
ometers along with their associated electronics. 
Outside the stable member are three servo-driven 
isolation gimbals (azimuth, elevation, and roll, 
in that order) isolating the instruments from air- 
craft rotation. The support structure surrounding 
the isolation gimbals is designed so that the 
tracker assembly mounts directly to its base, there- 
by providing a physical tie between the structure 
of the inertial reference system and the tracker 
pointing axes. 

The IMU has a unique thermal-control system which 
isolates the inertial instruments from the aircraft 
environment. The gyro mounting surfaces are held 
at approximately 46.1 ± 0.05°C and the accelerometer 
mounting surfaces are held at approximately 43.3 
± 0.05°C. This is accomplished by mounting the 
inertial instruments directly to the stable member 
air-cooled heat exchanger and employing thermal 
shims and individual controllers where necessary. 

An X, Y, 2 coordinate frame fixed to the inner, 
or stable, member of the IMU provides the attitude 
datum. This frame is called the indicated frame, 
platform frame, or "p" frame. It is initially estab- 
lished during a ground-alignment process, consist- 
ing of leveling and gyrocompassing . During leveling, 
the Xp - Yp plane is established normal to the local 
gravity vector. During gyrocompassing, the Yp axis 
is established normal to the earth's rotation vec- 
tor. The process of calibration and alignment of 
an inertial navigation system is a complex subject 
in its own right. A preliminary exposition of this 
process for the APT system is given in Reference 1. 

In the IMU, the sensitive axes of the accelero- 
meters are fixed in a known orientation with respect 
to the three coordinate axes of the "p" frame. . The 
outputs of the accelerometers are integrated twice 
in the computer to obtain position data for naviga- 
tion, and converted to latitude and longitude. The 
"p" frame is maintained when in flight by torquing 
signals proportional to the computed change in 
geodetic latitude and longitude. Thus, the plat- 
form is driven to maintain its attitude (except for 
errors) with respect to the normal gravity field of 
the earth and not the actual gravity field of the 
earth. It is this property of an IMU which allows 
the measuring of changes in the deflection of the 
vertical, a point discussed in detail in Reference 
2. In an error-free IMU, one would have a complete- 
ly self-contained instrument capable of giving posi- 
tion, velocity, and direction of the normal gravity 
vector in real time. 

The inertial system is a low-noise datum at high 
frequencies, capable of providing base-motion iso- 
lation from maneuvers or wind gusts. The inertial 


sensors (gyros and accelerometers) suffer from 
long-term drift, however, and cause errors in the 
low-frequency portion of the spectrum. To minimize 
these long-term drifts, the inertial components 
selected are high-performance instruments developed 
for an Air Force missile application. 

Even with the most accurate inertial systems, it 
is necessary to obtain updating from ground truth 
in order to reach the specified precision. To 
provide high-accuracy position updates, a laser 
tracker is employed. The tracker makes a vector 
measurement to a surveyed retroref lector ; range 
and two angles are measured. A minimum of three 
retroref lectors , not in a straight line, are re- 
quired in a 3- by 30-km area. The high-frequency 
information from the inertial system is optimally 
combined with the low-frequency information from 
the tracker to obtain an estimate much better than 
either. The statistical techniques employed for the 
optimum filtering of the data are given in the fol- 
lowing section. It can be shown that the important 
system errors, including position errors and plat- 
form misalignments, are independently observable. 

The laser tracker and laser profiler employ 
pulsed-laser rangefinding techniques. Distance is 
determined by measuring the round-trip time of a 
transmitted laser pulse. Pulsed gallium arsenide 
lasers are used as the transmitters for both de- 
vices. Constant fraction techniques are employed 
in the receiver threshold detectors to minimize 
timing errors due to received pulse-amplitude var- 
iations. In the profiler, pulse-selection tech- 
niques are employed to maximize the number of valid 
returns from terrain. The tracker contains a beam- 
splitter to separate the ranging and tracking func- 
tions. A four-quadrant detector is used to provide 
the error signals to the gimbal servos of the 
tracker. A precision time-interval counter is 
multiplexed between the profiler and tracker to 
provide the range measurements. 

A control and display unit enables the operator 
to sequence the system through its various operat- 
ing modes and to monitor and observe the status of 
the system. The system may also be sequenced auto- 
matically by the computer. Various system charac- 
teristics, including temperatures of inertial and 
other key components, critical voltages, and oper- 
ational modes are observed through the use of the 
control panel, as well as by computer monitoring. 
Steering information for the pilot to fly the 
desired courses for obtaining the required data 
is provided. 

Several military and commercial-type computers 
were investigated. For the APT system, the selected 
computer is the NORDEN PDP-11/70M, which has the 
speed, memory size, input/output capability, in- 
struction repertoire, architecutre , and software 
support necessary for the mission. 

Equations for Navigation and 
Topographic Profiling 

The first implementation of APT will be to prove 
the concept by performing only one mission, topo- 
graphic profiling, from the list given in the In- 
troduction. The same APT hardware, in the same 
aircraft, with the same on-board real-time software 
can perform the rest of the tasks by using differ- 
ent ground-based post-processing algorithms and 
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possibly different operational procedures. In this 
section the algorithms needed for the solution of 
the primary APT functions, navigation and profiling, 
are presented. 

When the aircraft reaches the ^survey area, two 
passes are made over three previously surveyed 
retroreflectors, deliberately chosen or placed not 
in a straight line. Range and angle data are ob- 
tained by the laser tracker in these passes by 
sightings on each retroref lector in sequence. These 
data are used to update and calibrate the position 
and orientation of the inertial system in prepara- 
tion for the profiling part of the mission. The 
positions of any unsurveyed retroreflectors pre- 
viously placed to satisfy any recognized need for 
added ground-truth references are determined rela- 
tive to the surveyed retroreflectors during the 
initial flight passes. Thereafter, each retrore- 
flector serves as a known reference point and may 
be tracked to obtain updating data as needed during 
the profiling phase of the mission. 

Figure 3 shows the flow of information through 
the onboard digital computer of the APT system. 

The input quantities shown are the profiler range; 
the three accelerometer velocity readouts and three 
gimbal angles from the IMU subsystem; and the two 
tracker gimbal angles and range to the retro- 
reflector from the tracker subsystem. The rates at 
which each of these quantities is sampled are in- 
dicated . 

The profiler subsystem data are sampled, time^- 
tagged, and subjected to reasonableness checks to 
eliminate signal dropouts and returns from the tops 
of foliage. The slope and intercept of the best- 
fit straight line are computed from 40 ms of those 
data which pass the reasonableness checks. The 
effect of the curve fit is to act as a low-pass 
filter of the profiler data. 

The IMU navigation subsystem samples accelero- 
meters for 80 ms (16 vector samples). The vectors 
are averaged and time-tagged at the mid-point of 
the data-acquisition interval. The effect of the 
averaging is to provide a simple low-pass filter on 
the input signals, in order to suppress further any 
aircraft vibration not taken out by the shock- 
mounting of the IMU assembly. Next, the accelera- 
tion vector is rotated by a fixed matrix to compute 
the north, east, and down components of accelera- 
tion, X, since the three orthogonally disposed 
accelerometers are mounted on the stable platform 
so as to receive equal components of gravity, g, as 
indicated in Figure 2. 

Equation (1) is the IMU navigation algorithm. 


IMU 


= A + g - 2W x R. 


IMU 


W x 


IMU 


( 1 ) 


where R is the position of the aircraft. 

This equation is solved, in the "p" frame above 

for R_,_. and R X __ I , by using trapezoidal 
IMU IMU 

£ule integration with a time step of 80 ms. The 
W vector is the angular velocity of the stable 
platform with respect to inertial space, given 
by Eq. (2) 


(fl + X) cos A 


W = 


+ A) sin A J 


( 2 ) 


where 0 = sidereal earth rate. 

A = the longitude rate. 

<j>,4> = the latitude and latitude rate. 

The g vector in Eq. (1) is a function of position 

g = gCRj^) and is given by the output of a gravity 

model; for example, the WGS72 ellipsoid. An error 
in position thus generates an error in g which, 
through integration of Eq. (1), generates more error 
in position. This feedback causes the characteris- 
tic 84-minute oscillations of error in the hori- 
zontal channels, and the characteristic doubling of 
the error in the vertical channel every 395 seconds. 
It is important to note, in this context, that one 
must use a gravity model to separate the earth's 
gravitational field from aircraft accelerations. 

The inertial instruments obey Einstein's fundamental 
principle of the equivalence of gravity and accel- 
eration in a reference coordinate system; they 
sense only the linear combination of the two. In 
what follows, we shall see that the laser tracker 
measurements, being referenced to the outside world, 
provide the mechanism for separating gravitational 
acceleration from nongravitational acceleration. 

The data from the tracker subsystem are sampled 
at high rates, low-pass filtered, and the tracker, 
navigation equation, Eq. (3)- is solved for rP . 
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In Eq, (3) the tracker range r is rotated about the 
y-axis for the inner tracker gimbal angle p; then 
rotated about the x-axis for the outer tracker 
gimbal angle 0; then rotated by a fixed orthogonal 
matrix C p for the alignment of the tracker assembly 
respect to the aircraft; then translated for the 
fixed distance, <f , from tracker assembly to IMU; 
and then rotated about the x, y, and z axes for the 
outer, middle and inner IMU gimbal angles, respec- 
tively. The result R P t is the distance from the 
retroref lector to the aircraft in platform (or "p" 
frame) coordinates. 

The IMU provides a slowly drifting (but smooth) 
navigation estimate at all times. The tracker 
provides a noisy but nondrifiting navigation esti- 
mate for short intervals of time (30 to 60 seconds) 
during which a retroref lector is within a 60° cone 
of visibility determined by the tracker viewing 
hole in the floor of the aircraft. As mentioned 
previously, the data from each of the navigators 
will be recorded for post-processing on the first 
flight trials. The algorithms in this section 
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which are used to process these data, however, are 
the same whether done in real time or after the 
fact. The rest of this paper proceeds with the 
description of the real-time implementation of data 
processing . 

The next step is to mix the data from the two 
subsystems. 


F = 


0 


F, 

/ 


0 


9 rows 


11 rows 


9 cols 11 cols 


where 


*> X -►p 

- R lt 


(4). 


By taking the difference between the positions 
indicated by the two navigators, as given by 
Eq. (4), the aircraft motion is subtracted out and 
there remains only a vector which contains a lin- 
ear combination of the errors in the two subsystems. 
The vectors D are then accumulated over a 2-second 
interval (25 vectors), and low-pass filtered again. 
The 2-second sample then represents input to a 
Kalman filter, or optimal estimation algorithm as 
'.escribed in Reference 3, Chapters 3 and 4. 

Kalman filter is specified very briefly by 
.lifting the elements of the state vector X, giving 
the system dynamics matrix F, the system measure- 
ment matrix H, and the measurement covariance 
matrix R, and the spectral density matrix Q. 

The state vector x consists of 20 elements. 

These are: 

■ ■+■ 

e = a linear combination of IMU position error 

and retroref lector survey error (3 elements) 

e = IMU velocity error (3 elements) 
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In order to build this F matrix, take W and W from 

Eq. (2) and R from the accelerometers. The notation 
used is that an underscore denotes the 3x3 anti- 
symmetric matrix associated with a given vector, and 


= 


B «= 
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vector of stable-platform misalignment 
angles (3 elements) 

a linear combination of accelerometer 
bias, vertical accelerometer scale-factor 
error, and gravity model bias error 
(3 elements) 

a linear combination of the constant com- 
ponent of gyro drift and the along-track 
gravity gradient errors (3 elements) 

a linear combination of various body-axis- 
fixed constant angle errors (3 elements) 

the bias error in the inner tracker gimbal 
(1 element) 


= 

R 


-1 


0 


0 


0-10 
0 0+2 


where R = radius of the earth 
e 

g = nominal magnitude of gravity 

In the' differential equation x = Fx + v, the (20 
element) vector v is called the plant noise. 
Spectral-density plots of laboratory measurements 
taken from actual APTS instruments revealed that a 
suitable model of the statistics of v is zero 
mean, Gaussian noise with a covariance matrix, 
called the Spectral Density Matrix, Q. 


6r 


the bias error in the tracker range 
measurement (1 element) 

Other elements which may be added to the state 
vector include: 

(1) The position coordinates of unsurveyed 
retroref lectors. 

(2) Unknown values of gravity disturbance vectors 
and gradients at a point where a retro- 
reflector is located. 

(3) The gravity gradients and even rate of change 
of gradients averaged over the flight traj- 
ectory. 

The present state of development of the software 
includes the above 20 elements, which suffice to 
account for all the APT instrumentation errors. 

This state vector x satisfies a linear differ- 

entxal equation x = Fx + v, where the system dynam- 
ics matrix F is a 20 x- 20 matrix which is parti- 
tioned as follows: 
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3 rows 
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It is possible to show (see Reference 4 page 64) 
that the vector D in Eq. (4) is, to first order, 
linear in the elements of the state: 
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i . e . , D = Hx + 

. . .where the time-varying H matrix has 3 rows 
and 20 columns. The vector n is the measurement 

noise in the laser tracker. Since five gimbal 
angle measurements and one range measurement are 
used to compute the tracker position (see Eq. (3)); 
since an angle error affects the tracker position 
perpendicular to the line from aircraft to retro; 
since a range error is along that line; and 
under the assumption that the angle errors are 
uniformly distributed; then, it is concluded that 

R = 5(n,. • n,J = a (uu ) + ■=• r a , (I - uu ) 
It It range 2 angle 


where a 


range 
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angle 


= the variance in range measurements^ 

> 

= the variance in angle measurements . 


= unit (R lt ) 
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This completes the brief description of the Kalman 
filter. 

Referrring again to Figure 3, the Kalman filter 

provides estimates x of the state and estimates P 

of the state covariance every 2 seconds. These 

estimates may be combine^ with the IMU data to 

produce a best estimate R , of aircraft position. 

a/ c 

If the aircraft has acquired data from at least 
three retroref lectors not in a straight line, and 


the quantities e, e, (the first nine elements of 
the state) are significantly different from zero 
(using a 3o test based on the square root of the 
diagonal of the covariance matrix), then Eq. (5), 
is used; otherwise, the best estimate is considered 
to be the uncorrected IMU data. 
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where e and e = state elements 

e = obtained from x = Fx 

t = the time of the most recent 

Kalman filter update, so that 

(t - t ) < 2 seconds 
o 

The corrected aircraft navigation data from Eq. (5) 
are used to build the state dynamics matrix F in 

the the Kalman filter, and W and W or Eq. (2). W 
is the gyro torque rates used to locally level the 
stable platform. 

Finally, the position of the aircraft R . from 

a/ c 

Eq. (5) can be combined with the profiler data to 
obtain elevation and position at the nadir. Since 
the profiler will be located several feet from 
the IMU, it will be measuring the elevation^of a 

point which is displaced horizontally from R . and 

a/ c 

which (in general) will thus have an elevation dif- 
ferent from the IMU, depending upon aircraft pitch 
and roll angles. The geometry of the aircraft 
gives (in locally level coordinates) 


^alt = K/c ~ ~ i) ’ c z (-“) ’ C y (-B) • C x (-y) • P b 

•y 

where R = position of the profiler. 

R# c = position of the IMU. 

C z (-a),C (-6),C (-y) are the same rotation 
v x matrices, dependent on 

IMU gimbal angles a, 6 
and y that are defined 
in Eq. (3). 

->-b 

p = the distance of the profiler from 
the IMU in aircraft body-axis 
coordinates . 

This leads to the actual survey datum: 


(x 


alt’ 


y alt’ Z alt + h) 


= R 


gr 


( 6 ) 


where x , .y ,^,z , are the north, east, and down 
alt alt alt ’ ’ 

components of $ 

Computer Simulations 


In order to verify the APT system concepts, there 
is an on-going effort to exhaustively model all the 
dynamics and statistics of a typical APT mission on 
a digital computer. The strategy is: (1) to write 

an off-line non-real-time version of the navigation 
and profiling estimation equations given in the pre- 
vious section, called an estimation model; (2) to 
write a more exhaustive, second-order simulation of 
the APT system error dynamics, called a truth model; 
and (3) to compare the outputs ^of the estimation 

■+ 

model (estimates of the state x) with the inputs x 
to the truth model, the residuals giving a realistic 
estimate of the eventual APT performance. Figures 
4 and 5 illustrate the use of the APT truth model 
and estimation model. At their current stage of 
development, both computer programs can give a 
realistic estimate of the amount of error contri- 
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buted by the APT instruments themselves. The pro- 
gramming for component of error contributed by 
residual gravity model errors is not yet complete, 
however, and so will be accounted for later in this 
section by another technique. 

In a typical computer run, we simulated a 20- 
minute section of a nominal APT mission with the 
truth model. The simulation flew the aircraft over 
the trajectory indicated in Figure 6, wherein the 
three circles represent the region of visibility 
around three perfectly surveyed retroref lectors . 

The retros are about 10 km apart, the aircraft has 
a mean speed of 53 1/3 m/s, a mean altitude of 

O. 91 km, and the viewing hole for the tracker de- 
fines a 60° cone of visibility. In this simulation 
the laser accuracy, after low-pass filtering for 

1 second, was 3 cm in range, 22 sec (about 0.107 
mrad) in pointing. The IMU performance was based 
on parameters and spectral models derived from 
classified military data taken on instruments iden- 
tical to the APT accelerometers and gyros. The 
truth model program simulated the aircraft dynamics, 
IMlI error drift, laser measurement errors, and 
retroreflector geometry, but not the laser profiler 
or the gravity disturbance vector. 

The 20-state Kalman filter described previously 
was used as an estimation model on the data gener- 
ated by the truth model. For each state the reside 
ual error was computed as a function of time, and 
plotted. The Kalman filter fluctuates wildly for 
the first part of the flight, until data from three 
different retroref lectors (not in a straight line) 
have been processed, and the filter achieves geo- 
metric sufficiency. Accordingly, Figures 7 and 8 
corresponding to the vertical position and acceler- 
ation error, respectively are plotted only for the 
last 700 seconds of the 20-minute flight. On each 
graph, three curves are plotted. A middle curve 
representing the state-element residual estimation 
error is enclosed by an envelope corresponding to 
±1.6 a (90% confidence), where a is computed by the 
Kalman filter as the square root of the appropriate 
diagonal element of the estimated covariance matrix 

P. 

In Figure 7, we plot the vertical position error 
in cm. Since the error curve fills, but remains 
inside, its 90% envelope most of the time, we may 
say that the Kalman filter is a self-consistent 
model when driven by this truth model data. Qual- 
itatively, the curve and envelope show character- 
istic periods of compression and expansion corres- 
ponding to the aircraft being inside or outside 
a retroreflector cone of visibility. When outside 
a cone of visibility, the envelope expands as time 
to the 3/2 power. 90% of the errors are less than 
8 cm. This figure represents the fundamental sur- 
veying error from the APT instruments; however, if 
surveyed retroref lectors were placed closer togeth- 
er, then we might expect even better performance. 

In Figure 8, the vertical acceleration error is 
plotted in mgals. It can be seen that the filter 
is still converging at the end of 20 minutes, how- 
ever, 90% of the data are already within a 0.1-mgal 
band. What we have shown here is the ability to 
determine the mean gravity anomaly over the region 
(actually averaged over the trajectory of the air- 
craft) to within 0.1 mgal absolute accuracy. More 
sophisticated modelling will allow determination of 
mean gravity gradients over the region, or mean 


gravity gradients between any pair of retro- 
reflectors or values of gravity disturbance vectors 
above each retro at the height of the aircraft, all 
to within an absolute accuracy of 0.1 mgal. Plac- 
ing more retroreflectors will extend the region of 
the gravity survey beyond the simulated 30 km; plac 
ing them closer together will increase the resolu- 
tion of a gravity survey better than the simulated 
1 point in 30 km. 

The other 18 state elements were plotted in a 
similar format, yielding information on navigation, 
horizontal retroreflector survey error, deflection 
of the vertical, gyro drift, APT system mechanical 
misalignments, IMU platform attitude error, and 
laser tracker errors. These results will be pub- 
lished in January 1979. 

So far only the APT instrumentation errors have 
been simulated. A back-of-the-envelope computation 
was carried out using gravity anomaly data from the 
Denver area supplied by the Department of the 
Interior, Geological Survey, to demonstrate the 
feasibility of correcting for the effects of grav- 
ity anomalies. The study consisted of calculating 
position errors attributable to residual gravity 
errors along a 30-km flight path by integrating 
twice. It was assumed that the aircraft would fly 
along the path at a constant velocity of 53.3 m/s, 
resulting in a 600-second flight. The initial 
position and velocity errors were assumed to be 
zero and these errors were again set to zero each 
time it was necessary to simulate the effect of 
three laser tracker updates. In addition, it was 
assumed that measurements of gravity anomalies were 
known a priori, or could be estimated, at three 
retroreflectors located along the flight path at 
13-km intervals. Based on these three measure- 
ments, a quadratic model for the gravity 'data was 
chosen and the coefficients of the model were com- 
puted. The difference between the known gravity 
anomalies and the gravity values obtained from the 
quadratic model were considered to be the residual 
gravity errors and thus the cause of position 
errors. 

The result of this study was that the standard 
deviation associated with the vertical position 
errors was 6.1 cm, the maximum error was 21 cm, 

90% of the vertical position errors were bounded 
by 12 cm. 

Although this study represents only a first cut 
at the gravity problem, it suggests that a combina- 
tion of modeling, estimation, and a priori measure- 
ments should be sufficient to account for gravity 
anomaly errors in the APT system. 

Since the vertical error from instrumentation 
is within a 90% confidence interval of 8 cm, since 
the unsimulated gravity disturbances will be about 
12 cm, and since these are uncorrelated, then 

RSS error bound = 4 8 2 + 12 2 < 15 cm 
Crustal Movement Monitoring 

In the preceding section, the capability of the 
APT system to perform real-time navigation and 
topographic profiling to an accuracy of 15 cm in the 
vertical, relative to three control points, was 
demonstrated in a survey area about 30 x 3 km. For 
those applications where the survey region has been 
chosen to straddle a fault line and where one plate 
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has moved relative to another by more than 15 cm, 
the APT system could be used, unmodified, by the 
simple expedient of placing all three control points 
on the same plate; differencing the profiler eleva- 
tions taken on an APT survey at time tg from an APT 
survey taken at time t^; and averaging the differences 
on each plate separately to remove bias errors and 
the effect of not repeatin g the aircraft path ex- 
actly. If Zj(t k ) - 2 j ( t q ) is the average of APT 
vertical displacement measurements of plate J at 
time t k , then the vertical movement of plate 1 with 
respect to plate 0 during the time from tg to t^ is: 

[ Z i(ti> - Z^tg)] - - Zg(tg)] 

There are some problems with this method: resolu- 

tion, extent of survey area, and accuracy. By 
modifying the APT system slightly we can solve these 
problems . 

First of all, the resolution problems can be 
solved by eliminating the profiler. Profiling data 
is taken every 40 ms (i.e., 3 meters) then has to be 
averaged over the aircraft trajectory to give a 
resolution meaningful for crustal plate dimensions. 
Also, profiling data will not be exactly the same 
in successive surveys, since the aircraft can not 
repeat its flight path exactly; hence, the survey 
area has to be profiled densely so that the averages 
Zg(t^) are representative of the area. Lastly, the 
laser tracker can be used as a more effective sur- 
veying device than the profiler. This can be accom- 
plished by inverting Eq. (3) and (4) and solving for 
R rr . Instead of using measured tracker range and 
pointing angles with a known retroreflector position 
to compute the aircraft position, the aircraft posi- 
tion is assumed known and the rotation and trans- 
lations performed to compute the location of the 
retroreflector . The accuracy will be as good as 
the IMU subsystem navigation accuracy. The air- 
craft needs to repeat its trajectory from survey to 
survey only well enough to enter the cone of vis- 
ibility of the retroref lectors . Retroref lectors 
can be placed atop fixed monuments in a grid at 
whatever spacing corresponds to meaningful resolu- 
tion of crustal movement measurements. If Rj(t^) 
is the APT surveyed position of retroreflector 
J at time t k , then the movement of retro J with 
respect to retro I from time tg to time t^ is 

[Rj( tl ) - Rj.^)] - [RjCtg) - Rj.(tg)] (7) 

Notice that we compute vector displacement, not 
just the vertical component. 

The problem of survey extent was related to sur- 
vey accuracy. As mentioned in the previous section 
in reference to Figure 7, position errors in the 
IMU grow as time to the 3/2 power when there are no 
tracker measurements. One way to increase the size 
of the survey area is to fly faster so that more 
distance can be covered for a given error bound. 

In a jet aircraft, speeds of 266.7 m/s, or five times 
the speed assumed in the previous section, may be 
assumed. Such speeds imply a higher safe operating 
altitude, say 4.5 km. 

The problem of accuracy was addressed by a com- 
puter simulation of the APT system in which an air- 


craft was flown at a height of 4.5 km, at a speed 
of 266.7 m/s, repeatedly around a closed, oval tra- 
verse about 120 km long and 30 km across. At each 
end of the oval was a cluster of three surveyed 
retroreflectors about 20 km apart in a right tri- 
angle configuration. The clusters were about 100 km 
apart, and the complete traverse could be flown in 
20 minutes. The parameters for the laser tracker 
were the same as for the profiling mission. If the 
navigation error from this simulation is observed, 
we get a measure of the ability of APTS to survey 
in retroreflectors. In this simulation we did not 
include the effects of unmodelled gravity errors on 
the system, because these may be reduced to resid- 
uals by availability of a priori gravity surveys, 
reduced further by the ability of APTS to estimate 
the gravity field, and removed through the following 
survey operation. If we fly the same (±500 m) tra- 
jectory at the same speed, then the gravity- induced 
IMU errors will be the same from survey to survey 
and subtract out in Eq. (7) . Of course, the local 
gravity field may undergo temporal change also; 
however, the effect, which is 3.3 microgals per cm 
of vertical displacement is distinctly second order. 

The result of the simulation was that 90% of the 
errors were less than 18 cm. Such a result shows 
that the APT could be used for measuring only very 
large vertical crustal displacements unless further 
modifications are made. One improvement would be 
in more sophisticated modelling to use all of the 
information available. There is a post-processing 
technique, called back smoothing, (see Reference 3), 
which is a refinement to the Kalman filter algorithm 
offering the potential of significant improvement. 

A forward Kalman filter forms the best estimate of 
the state vector given initial conditions, the 
models, and all the measurements up to the time of 
the estimate. A backward Kalman filter forms the 
best estimate of the state given final conditions, 
the models, and all the measurements which occurred 
after the time of the estimate. An optimally 
weighted average of the forward and backward filter 
outputs then forms the optimally smoothed estimate 
of the state. Under the assumption of error 

growth, it is easy to show that the peak error the 
optimally smoothed estimate is 1/4 the peak error 
in the Kalman filter estimate. The back smoothing 
algorithm is very cumbersome, costly of computer and 
programming time, and yet offers the possibility of 
measuring crustal movement with the APTS to within 
about 18 t 4 = 4.5 cm. 

We shall now discuss modifications to the APT 
system hardware and improvements to mathematical 
modelling which would push the APT performance to 
its ultimate limit. What has been shown, so far, 
is that the errors derived only from APT instrumen- 
tation are about '2.5 times higher for measuring 
crustal movements (_>_ 18 cm before backsmoothing) 
than for topographic profiling. This factor occurs 
because the geometry of the tracker measurements 
is so much less favorable at the higher altitude 
and sparser retroreflector spacing used to mea- 
sure crustal movements. This geometry affects the 
system through the t^/2 error growth between 
retrosightings and through the degradation in 
tracker accuracy at higher altitude. The time 
between tracker position fixes could be shortened 
only by flying even faster, involving the choice 
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of a supersonic aircraft; however, it was felt that 
aircraft maintainance expenses for such a vehicle 
were not. justified by the marginal gain in APT per- 
formance. 

The degradation in tracker performance with 
increased altitude occurs for two reasons. First 
of all, the laser beam must pass through more of 
the atmosphere, and over a greater range of pres- 
sure, temperature, and humidity. Therefore the indej 
of refraction will vary over a greater range along 
the ray-path, hence (through Snell's Law) there will 
be more atmospheric ray bending. Because of ray 
bending both the direction and distance measurements 
of the tracker are in error. If the APT system were 
modified so that the tracker were replaced by a 
two-color laser, then one could correct for the 
change in arc length due to raybending. The scheme 
is described in Reference 5. The pointing uncer- 
tainties can be reduced to 3 sec by atmospheric 
modelling of 60 (raybending distortion) versus 0 
(angle of incidence) as described in Reference 4. 

The second effect is a purely geometric one. The 
distance error caused by tracker pointing errors 
is 1 mm per microradian for each km of altitude. 

The precision angle readouts, as envisioned for 
APTS, have a = 22 sec effective pointing error. At 
0.91 km altitude this corresponds to 10 cm, at 
4.5 km to about 0.5 meter. Most of this error is 
repeatable and could be removed by careful labora- 
tory calibration of a table of 60 versus 0 for each 
precision angle resolver. If this were done the 
effective pointing error could be reduced by at 
least a factor of 4 to 5.5 sec. 

In the computer simulation which estimated 18-cm 
accuracy for Eq. (7) (4.5 cm presumed after back 
smoothing) we assumed all tracker information avail- 
able for the 100-km flight path (37.5 seconds fly-r- 
ing time) between the two clusters of surveyed 
retros was used for survey and none for navigation. 

In fact, such navigation information is available, 
and could be used if more sophisticated modelling 
were performed; thus, the results of the computer 
simulation are unduly pessimistic. There are three 
major sources of navigation data not included in 
the simulation. First, each retro supplies a vel- 
ocity measurement with which to update the drifting 
IMU. The 1600-Hz laser tracker pulse rate gives a 
complete time history of the aircraft with respect 
to, the retro, and that retro is not moving for the 
10 seconds of the survey. In fact the measurement 
equation for a velocity update can be obtained by 
differentiating Eq. (3) and (4). If one places 
enough retroref lectors along the survey trajectory, 
then one simultaneously surveys-in more points for 
measuring crustal plate movements and supplies 
almost continuous velocity information to the system. 
The limit is reached when the retros are so close 
that the regions of visibility overlap, because it 
takes 2 or 3 seconds to move the tracker from lock- 
ing on one retroref lector to the next. The effect 
of such a near-continuous velocity update on posi- 
tion accuracy should be quite dramatic, since, in- 
stead of a characteristic t3/2 error growth caused 
by twice integrating random acceleration errors, one 
would expect tl/2 error growth caused by a single 
integration of the nondrifting random residual 
velocity errors. In Reference 4, Chapter 3, Appen- 
dix B an estimate of less than 1 cm for t = 375 
seconds is made assuming a perfect and uninter- 


rupted velocity update. 

The second major source of navigation information 
comes from the fact that the crustal motions of 
closely spaced points will be highly correlated. 

The correlations, if modelled, would allow the 
extraction of position information for navigation 
from redundantly placed retroref lectors . 

The third source of information occurs because 
there is a bandwidth separation of the navigation 
error and APT survey error which has not been taken 
advantage of. So far, we have been observing navi- 
gation error as though this were the same as survey 
error. However, the navigation errors are pure 
functions of time, and the survey displacements 
which we are attempting to compute are pure func- 
tions of latitude and longitude. If one were to 
repeatedly cross over a retroref lector , then one 
would expect the IMU error to be statistically 
independent from crossing to crossing, so that the 
survey error could be a factor of /n times better 
than the navigation error, where N is the number 
of times the retro has been crossed. 

Unfortunately, updating the existing APT simu- 
lation program with the equations necessary to 
implement velocity updates and survey computations, 
does not fit the APT development schedule, prior- 
ities, timing, and budget. As mentioned in the 
previous section, the system will be programmed to 
perform the tasks of navigation and profiling. What 
can be done, however, is to use the existing simu- 
lation program to infer a best case bound on APT 
performance. There is already a worst case bound, 
4.5 cm with back smoothing, gotten by assuming no 
tracker information available during an APT survey. 
The best case bound is gotten by using the existing 
computer simulation to supply position data from 
perfectly surveyed retros and observing navigation 
error. Admittedly the combined effect of velocity 
updates, survey correlations, bandwidth separation 
of survey and navigation errors, and back smoothing 
can not be as good as position updates, but it does 
set an APT performance limit. 

Accordingly, another computer simulation was run, 
with the IMU model the same as APTS, with the laser 
tracker pointing accuracy improved to a = 5.5 sec. 
Sixteen perfectly surveyed retroref lectors about 
16 km apart with a 60° cone of visibility lay along 
an oval flight trajectory as indicated in Figure 9. 
The aircraft was flown at 266.7 m/s at an altitude 
of 4.5 km. The vertical position errors for one 
traverse of the oval are plotted in Figure 10 in 
the same format used in Figure 7. In Figure 10, 
one can see by inspection that most of the errors . 
lie between ±1.8 cm. This number was confirmed by 
plotting the points in Figure 10 in an integrated 
histogram format. In Figure 11, the abscissa is 
error bound in cm, the ordinate is the % of errors 
from the simulation less than or equal to the 
abscissa. From Figure 11 one can read off 1.8 cm 
as the 90% error bound. 

From these computer simulations we have reached 
the conclusion that the APT system could survey 
vertical crustal movements much better than 4.5 cm 
relative to three control points, but cannot sur- 
vey better than 1.8 cm, in a survey area of di- 
ameter 120 km (about 1.5 - 4 parts in 10 . The 
system also has the capacity to measure horizontal 
displacements to the same order of accuracy. 
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Figure 1. USGS/APTS functional diagram. 


I 



Figure 2. APTS Inertial Measurement Unit (IMU) and 
laser tracker schematic. 
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Figure 8. Vertical acceleration error and covar- 
iance envelope from navigation simulation. 
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Figure 10. Vertical position error and covariance 
envelope from crustal-measurement simulation. 



Figure 9. Flight trajectory and retroref lector Figure 11. Integrated histogram of vertical posi- 
geometry for simulation of APTS measurement of tion errors from crustal movement simulation, 

crustal movements. 
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